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This paper investigates the density expansion of the thermodynamic properties of a two com- 
ponent plasma under the influence of a weak constant uniform magnetic field. We start with the 
fugacity expansion for the Helmholtz free energy. The leading terms with respect to the density 
are calculated by a perturbation expansion with respect to the magnetic field. We find a new mag- 
netic virial function for a low density plasma which is exact in quadratic order with respect to the 
magnetic field. Using these results we compute the magnetization and the magnetic susceptibility. 
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I. INTRODUCTION 

The topic of matter in magnetic fields has received much attention. The magnetic field has a great effect on the 
individual motion of charged particles. In classical mechanics the motion of the free particles in a magnetic field can 
be described by circular orbits. The frequency associated with the rotation is the cyclotron frequency Uc — eBo/nie. 
In quantum mechanics the motion perpendicular to the magnetic field is quantized with the corresponding energy 
eigenvalues E± = hcjcin + 1/2) A wide range of subsequent investigations covers the properties of atoms and 
molecules in magnetic fields. Their study is motivated by the astrophysical implications concerned with the physics of 
pulsars and neutron stars, but also by its application in quantum chaos [^Q. The calculation of the energy spectrum 
of a hydrogen atom in strong magnetic fields has been tackled by various authors One of the complications 

found in the theory is the coupling of the center of mass motion with the relative motion |^ . This also complicates 
the calculation of the thermodynamic functions of a low density plasma. However, this effect becomes important 
at strong fields only. Throughout this paper we will consider the weak-field limit in which case this effect becomes 
negligible. In this limit the proton mass is considered to be infinite. 

Although the magnetic field affects the individual motion of the particles, there is no influence of the magnetic field 
on the equilibrium properties of a classical charged particle system. This follows from the Bohr-van-Leeuwen theorem 
1^,^, which can be easily derived by changing the variable in the momentum integrals in such a way that one works 
with the variable tt = p — eA. As a result of this all equilibrium properties are independent of the magnetic field. In 
quantum mechanics this argument is no longer valid, since the momentum operator and the coordinate operator of 
a particle do not commute. A common example of an equilibrium value which depends on the magnetic field is the 
magnetization of an electron gas. It was shown by Landau pO[ | that the low field magnetization of a spinless electron 
gas, in Boltzmann statistics, is 

M... ^ , (1) 

12m 

i.e., the so called orbital part contributes to a diamagnetic response. However, the full response of the free (noninter- 
acting) electron gas, including the spin part, is paramagnetic. 

The magnetization of a system of charged particles is a boundary effect. In classical mechanics the magnetization 
induced by the motion of the bulk electrons is cancelled by the magnetization connected with the surface current. 
Again in quantum mechanics this statement is not valid anymore. Landau used a perturbation expansion of the 
free energy with respect to the magnetic field, to circumvent difficulties due to the boundary effects. In doing so 
the electrons at the boundary of the system were neglected. Then the magnetization is found as the derivative of 
the free energy with respect to the magnetic field. Another approach has been chosen by Teller |jll[. He calculated 
the current at the boundary of the system, produced by the motion of the elctrons under the influence of uniform 
magnetic field. From this he computed the magnetization of the system and was able to show the equivalence of his 
method and Landau's approach. However, as Teller already pointed out Landau's method is much better suited for 
more complicated problems. 

The difficulties connected with the boundary effects are perhaps one of the reasons for the few results concerning the 
equilibrium statistical mechanics of a low density quantum plasmas embedded in an external magnetic field. A first 
attempt beyond Landau has been pursued by Alastuey and Jancovici. They studied, by means of a Wigner-Kirkwood 
expansion, the magnetic properties of a nearly classical one component plasma (OCP) in two and three dimensions 
in the weak field [Oj as well as in the strong field limit [Ol . Related problems were treated by Cornu 13 and Boose 
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and Perez ]15| who derived a formally exact virial expansion of the EOS of a multicomponent plasma by using the 
Feynman Kac path-integral representation of the grand-canonical ensemble. 

This paper is aimed to calculate the magnetic properties of a quantum plasma in the low density limit. These 
systems are characterized by a small coupling parameter F, which is given by 

^ " AneokTd ' 

where d = (3/47rn)^/'^ is the mean distance between the particles. We follow the method of Landau for the calculation 
of the magnetization. The starting point is the fugacity expansion of the Helmholtz free energy. In a previous work 
p6[ the authors have performed a perturbation expansion for the equation of state of a low density plasma up to 
the order n^e^, which is valid at arbitrary magnetic field strength. In the present paper we perform an expansion 
of the thermodynamic functions for a low density plasma (F < 1) up to the order n? and calculate the coefficients 
of this expansion in quadratic order with respect to the magnetic field, without making any approximation with 
respect to the Coulomb problem. Using a diagrammatic language this can be restated as the calculation of all ladder 
diagrams expanded to second order in the magnetic field. In doing so we will consider the B and terms of the 
Hamiltonion separately. We will see that the separate contributions are divergent, only the sum of all contributions 
give a convergent expression. This is the price which we have to pay within the present method circumventing the 
calculation of boundary effects. 

This paper is organized as follows. In section II, we discuss the representation of thermodynamic functions by a 
fugacity expansion. The second virial coefficient of a magnetized plasma will be discussed in section III. In the first 
part of section III, we present in more detail the calculation of the electron-ion contribution to the thermodynamic 
functions in the case of an infinite proton mass and in the second part an analytical continuation will lead us to the 
electron-electron contribution. In the third part of section III the asymptotic behavior of the new proposed magnetic 
virial functions will be studied. Finally, the derived results are used to compute the magnetization and the magnetic 
susceptibility in section IV. 



II. REPRESENTATION OF THE THERMODYNAMIC FUNCTIONS BY A FUGACITY EXPANSION 

In this section, we briefly present the general method used in this work and give the exact results derived in an 
earlier work . Let us consider a two-component charge-symmetrical system of N spin half particles of charge (-e) 
and mass rrie and N spin half particles of charge e and mass rrii . The Hamiltonoperator of our system consists of two 
particle contributions. Each pair of species a and b contributes 

Hab = + /^^oCTz + ^ + ^J■BBocrz + A\4h(r) , 



with the Coulombic interaction potential 



= -!,+! (3) 



K.(r) = . (4) 

Here Hab is the Hamiltonoperator of the two particle system and iJ^^ of the noninteracting system. The additive term 
fi'^Boaz takes into account the coupling between the intrinsic magnetic moment (/i^ — Cah/ {2ma)) of the charged 
particles and the magnetic field. 

We suppose that the pressure can be split into ideal and interaction contributions 

P^Pid+ Pint ■ (5) 

In the case without Coulomb interaction = the pressure and the particle density of the plasma in a homogeneous 
magnetic field B = (0, 0, Bq) are given by a sum of Fermi integrals over all Landau levels n 

a ° Ti=0 a °- n=0 



{xa = huj^/{2kT) with w° = \ea\Bo/ma, Aa = h/ ^j2TTmakT , and z° = exp — nhuol)]). The prime indicates the 
double summation due to the spin degeneracy except for the n — Q level. 



2 



The interaction part of the pressure will be expressed in terms of a fugacity expansion which will be truncated after 
the second virial coefficient ||l7H2l[] 

^^P^nt = ^+Y.~^ah (^Xl.eabH^^ab) + ^^J^^ + Ba^ + Qiz"- \n ~z) , (7) 



26 

where we have introduced the modified fugacities 



° tanh(a;a) ' 



(8) 



in order to have Za ria in the limit of small densities. The first term on the r.h.s of Eq.(|7l) is the Debye contribution 
in the grand canonical ensemble. The squared inverse Debye radius is given by = /3(£r/eo)(2e + Zi). Since it is a 
classical contribution the Debye term does not depend on the field. In the limit of small densities Eq.(0) coincides 
with the formally exact virial equation of state derived by Cornu |Q and Boose and Perez . We now try to extend 
these calculations and focus on the calculation of the second virial coefficient Bab- In order to avoid convergence 
problems let us in a first approach cut the Coulomb tail at large distances, i.e. Vabir) = if r > i?. Then the second 
virial coefficient reads 

^ ^ / A|tanh(xJ\ / A| tanhM \ ^ _ ^_0o^^ 
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This function will be studied in more detail in the next section 



III. SECOND VIRIAL COEFFICIENT OF A MAGNETIZED PLASMA 

It is convenient to divide the second virial coefficient into the direct part _B^^ and the exchange part B^a 

Ba, = B^ab + BTaSab, (10) 

and to compute them separately. The exchange part of Bab as defined by Eqs.(^) with the Hamiltonoperator given 
by Eq.(^) is convergent while if i? — > cxd the direct part is divergent. This is due to the long range behavior of the 
Coulomb interaction which leads to collective effects. In order to include these collective effects one has to perform a 
screening procedure, which may then lead to convergent expression for Bab- Such a technique is well established in 
the zero magnetic field case p7|-pl| and can be easily extended to the nonzero magnetic field case. In general, Bab is 
an analytic function of the interaction parameter ^ab pl| , defined by 

Cafc = ~"; TTfTT 1 (11) 
AttcqkI Xab 



with Xab = h/\/2mabkT and mab = Tnami,/ {ma + mb) beeing the effective mass. Hence Bab may be written as a 
Taylor expansion. Using the methods as described in Jl?]-^ we derived in our earlier work the lowest order 
contribution to Bab- As in the zero field case we write the direct part of the second virial coefficient of the plasma in 
the following form 

Bab = B'ab + B'ab , (12) 

where the contributions of second and third order in ^ab arc included in -B^j, with 



B'ab = -l^^^^>^lbCbh2{Xa,Xb) " ^ ^ ^'Ib^lb^siXa, Xb) 



(13) 



In general the magnetic field correction ft,2.3 satisfies /i2,3 = 1 if the magnetic field B = Q. The second order term has 
been found in p6[ and is explicitely given by 



h.ix.,Xb) -iU-fdt V^) iVa + Vb) ) , (14) 
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with ya,b = '^^aabb svah{xa,bt) sinh(a;a,b(l — 0)/(-^a6*(l ~ t)'^Xa,b sinh(xa,fc))- The magnetic field correction /13 is so 
far not exactly known. In the limit of zero field /13 = 1 holds; furthermore, in the next section we will derive an 
expression for /13 in the weak field limit. 

Formally the higher order contributions may be expressed by a resolvent expansion pO[ | 

I f tanh(a;a)\ / tanh(a;h)\ 



2n[2 xa j I 2 Xb ) ^" 



TrY—f dze-"' 



^ Vab 



ab 



' 1 

ab 



The contour integral may be taken in the sense of an inverse Laplace transform. The operator P" means that all 
terms of order less than have to be omitted, since they have already been taken into account in B'^i^. The series 
may then be written in the general form 

0^3/2x3 Cjk ~ '2,)hkiXa,Xb) f^ab\ 

Bab-^n^ x^,}_^^ r(i + /fc/2) [-YJ ■ ^^^^ 

The functions hk expressing the magnetic corrections satisfy the zero field condition 

/ifc(0,0) = l. (17) 

Therefore in the zero magnetic field case an exact calculation of the convergent second virial coefficient is possible in 
agreement with earlier work [20|pl[| . 

An alternative expression for the field free virial coefficient which we may refer to as SjJ^ may be obtained by 
introducing the quantum virial function Q^{£,ab) according to 

S^, = 27rXl,Q^iU) , (18) 

with 



Q°iU) = -Iv^eab - Ulb + - + rTT 

^ ^ fc=4 ^ 



C(fc-2) fU^" 



k/2) 



(19) 



Note that the second order term may be included into the series, since C(0) = ^1/2. 

Now let us discuss the exchange part. Again, as it was shown in Ref. [[l6| this contribution may be written in a 
Taylor expansion 

Baa=-' Ka^{^-'i ) r(l + ^2) \ — ) ' ^^^^ 

Here we have included the terms with nonpositive arguments in the C function. In particular we have used the relation 

lim^(l - 22-'=)C(fc - 1) = ln2. (21) 

The zero field results are reproduced, since we have 6^(0) = 1 and they may be written, after introducing the exchange 
virial function i?°(^aa), as 

B^T^-^Xl^E^U). (22) 

with 

The influence of the magnetic field on the exchange part has been studied in for the lower order terms 5o, &i, and 
&2 • The following analytical expression were derived 
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cosh(2a;a) tanh(a:;a) 
Oo[Xa) — , 

cosh (Xa) Xa 



and 



arctanh-v/ 1 ^ 



^ ^ _ cosh(2a;a) tanh(a:a) y ^„ ^25) 

COsh^(Xa) Xa ^/^ tanh(3:„) 



For an integral representation of 62 we refer the reader to Eq.(C4) of p^ . 



A. Expansion in the weak field limit for the ion-electron interaction 

We shall find an expansion of the second virial coefficient in the weak-field limit. The magnetic field is now assumed 
to be a small perturbation to the field free Coulomb problem. In this case we can use the already established results 
for the second virial coefficient in the absence of a field Due to the invariance of the thermodynamic functions 
under the transformation B — > — B the first correction to the field-free results will be quadratic in the magnetic field. 
This can also be verified in the ideal contribution (Pauli spin magnetism and Landau diamagnetism). 
Let us consider a hydrogen plasma with an infinite proton mass. This is a reasonable approximation in the weak 
field limit, as the proton frequency is proportional to the inverse mass of the proton. We chose the symmetric gauge 
A = 1/2 (B X r). Then the Hamiltonian in relativ and center of mass coordinates takes the form 

2mi 2me 8 2 ocp ^ire^r 

where is the electron cyclotron frequency. The elctron-ion contribution to the second virial coefficient is given 
by the following trace 

= i^'" (f ) {f-^) ^ - ^'"^'^ (2^) 

As in the zero magnetic field case we have defined an operator P''" that takes into account the divergency, by omitting 
all terms of order with k < k' . 

The trace over the center of mass coodinates and over the spin variable is readily carried out. Again we use the 
resolvent representation to obtain the following contribution 

with Ae — hj \j2mfkT . Here we have introduced the Hamiltonoperators for the free relative motion 

and for the relative motion of the interacting particles 

h.. = hl-^^. (30) 

We are interested in the case of a weak magnetic field without making any approximation with respect to the Coulomb 
problem. For that we expand B^i in powers of x^. = TujJcl2kT. It can be easily shown that the linear term is equal to 
zero and the first nonvanishing term is proportional to x^. This contribution may be written as 

.J..(,4)..,..,..V (31. 

The first term comes from the expansion of the normalising constant and therefore the trace is solely given by the zero 
field result . In order to take into account the infinite proton mass in i3°j one has to replace Aei by Ae = ?i/ \/2mekT 
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and, hence, ^ei by = — eeei/(47reofcTAe). The other two terms, beeing of the order 0{B'^), are the result of an 
expansion of the trace in powers of the magnetic field and read as 

Bl, = 4.3/2;^3p.' / ( (-^^il) T^] (32) 

Jc^TTi \h''-z \ 2 dct)J h'^-z \ 2 d(j}J h^-zj ^ ' 

and 

Bi ^ -4,V'A»P'' I i|e-IV ^) , (33) 

Here h'^ = p^/2me — e'^/Aneor is the Hamiltonian for the Coulomb problem for zero magnetic field. In what follows we 
briefly outline the steps leading to the final result for Bl- and B"^- . For simplicity the calculations of these contributions 
may be carried out separately, but as will be seen below only the sum of both gives a convergent contribution. 

1. Calculation of 

Let us first concentrate on the calculation of B^^. The perturbation operator has spherical symmetry. Thus it 
is convenient to use the eigenf unctions of the Coulomb operator. With that the calculation of the matrix elements 
becomes trivial. As in the zero field case we can write 



m~ 



I n=l ;=0 m=-/ V-^" ^1 



f 4 



1=0 m=-l 



Here we have made use of the relation between the density of states for the continuum states and the scattering 
phase shifts 5i{k) of the Coulomb system. The eigenvalues of the Coulomb system read as En = — l/2n^ and can be 
expressed in terms of the parameter ^ by 

(35) 

First we compute the discrete part of the partition function, that is given by the first term in ( [3^ ) and reads 

2 oo n — 1 / 

^ ' n=l 1=0 m=-l 

where we have performed the inverse Laplace transform. The summation over m and 1 is trivial and one immediately 
finds, that 



Bl^ = ^7:^^X1^ F^' Y: {n^ " n') exp . (37) 



12 ^ ' \2n 

n— 1 ^ 

By expanding the exponential and using the representation of the ^-function we obtain 

_ 

2 



Rlfc ,^V2^^^lr.k' C(fc-4)-C(fc-2) (ie 

Se. =4.^A,-P ^ r(l + k/2) U 



(38) 



So far we have calculated the bound state contribution to B],^. In the next step we consider the contribution of 
continous spectrum. For that we need the scattering phase shifts of the field free Coulomb problem that are given by 



d . . If e nie 



(39) 
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Making use of this relation and introducing t = X^k the second term in Eq.(p4) may be written as 



In order to compute the sum over m and 1 we will expand the ?/;-function, we have 

i?e^fz + l + -^|ee|) =i?ef;^V('=)(/ + l) f-^ICel)' . (41) 



Now the summation over m and 1 may be carried out. We obtain by introducing r = Z + 1 + s 

OO r — 1 

= (-l)^P'='i(C(fc~3)~C(fc-l))- (42) 
Next we perform all remaining integrations. In this context we may use the following integral representation 

c2m^ Jc,t^n^t-zf\2Vt^^'^) r((fc + 3)/2) 1, 2 ; ' ^^'^ 

Here the contour integral C in the complex t-plane encircles the positive real axis in the mathematical positive sense. 
Using Eqs.(|0|j4|andg|) we obtain, after shifting the summation index fc ^ A: — 1, the series 



Rl. o_3/2,3^e V^ C(fc-4)-C(fc"2) / ICelV , ... 

Finally we sum up the bound state (^) and the scattering state ( ^ ) contribution, which gives 

^i^2.3/^,3:^yC(fc-4)-C(fc-2) /u^y 

r(l + A:/2) V 2 y ^ ^ 

Here the sum runs from fc = 6, since in this derivation the lower order terms fc < 6 would give divergent contributions . 
However, formally the C function can be extended to negative values and therefore the sum to smaller k values such as 
fc = 2, 3, 4, and 5. It will be shown below that this extension is possible and gives the exact lower order contributions. 

Note that the bound state contribution and the scattering state contribution differ by a factor of 2. This general 
statement has been previously derived in the zero field case It is essentially a consequence of the analyticity 

of the second virial coefhcient BabiO ^^^^ expresses the fact of compensation of bound state and scattering state 
contributions according to Levinsons Theorem pC| ]. One may also regard it as rule of obtaining scattering quantities 
from bound state quantities. We will employ this relation in the following section. 

2. Calculation of B^^ 

Again we first concentrate on the calculation of the bound state contribution. Wc may use the eigenfunction of the 
Coulomb operator to evaluate the trace. Thus we have 



OO n — 1 I 



J /I J. I ^ 



nl 



m 



nlm ) . (46) 



The calculation of the matrix elements is readily carried out |i| , with the result 
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nlr 



m,u}, 



c 2 -2, 



ilr 



"^e^c 2^ (r 2 , . OKI ,\\ 

— ^aBy(5n +1-3Z(/-1)) 
Al^ + + 2{2l + l)m2 ~2l~2 



{21 + 1){21 - 1){21 + i) 

With that we obtain after integration and summation over the magnetic quantum number m 

n—l l—O 

By summing over 1 we get 



(47) 



(48) 



±:1 6Ci ' 



(49) 



As before we expand the exponential, introduce the (^-function and obtain the following expression for the bound state 
contribution 



x| ^ 7C(fc - 4) + 5C(fc - 2) (i, 
24 



fe=6 



r(2 + fc/2) 



(50) 



Now we shall calculate the scattering part. This contribution may be obtained by applying the same arguments that 
have led to the final expression of B\^^. 



o2 o_3/2 ,3^e V- 7C(fc-4)+5C(fc-2) / V 

i^e^- ^24^ r(2 + fc/2) \2) 



(51) 



Again, only contributions fc > 6 are retained from this sum. 



3. Final results for the electron-ion contribution 

We can now take the sum of the various contributions Eqs.(^, |4^ and pT] ) in order to obtain the quantum virial 
function. As we have indicated before we may now drop the operator P*"' and may postulate the virial coefficient, 
with 



Bi = 2nXlQ\i,) + 2^A^||g^(ee) , (52) 



,2 



where we have defined the new magnetic quantum virial function by 

-2) + (fc- 
r(2 + fc/2) 



^ (fc - 3)C(fc - 2) + (fc - 5)C(fc - 4) (UV 



k=2 

In spite of the fact that the derivation given above is valid only for fc > 6 we have extended the sum to fc > 2. By 
studying the asymptotic properties of this function we will show that the magnetic quantum virial function has the 
correct asymptotics for large ^. Another independent verification of this result can be obtained by expanding the exact 
second order contribution as given in Eq.(^). The quantum virial function Q^{^e) may be interpreted as the limit of 
Q^(Cei) with an infinite proton mass to^ — s- cx) . For fc = 3, 5 we make use of the relation limfc^3(fc — 3)C(fc — 2) = 1. 

In the next section we will show that the same analytical function determines also the contribution of the electron- 
electron interaction. 
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B. Electron-electron contribution 



We first study the Hamilton operator in cm. and relative coordinates. In this case the hamiltonian is separable 
and may be written as H^e = H"J^ + h^e , with the center of mass hamiltonian 

^eT = :^ + ^i?^sine^-.^|-. (54) 
Ante Anie 2 9$ 

It describes the free motion of a particle with mass 2me parallel to the field. While we have an harmonic oscillator 
with frequency uj^ = eB/m^ perpendicular to the field. The hamiltonian for the relative motion is given by 

e^B'^ o . ^, .fiujc d f? 

/lee = — + -7^r2sm02_j^ +^ . 55 

It has the same structure as the hamiltonian for the relative motion of an electron in the field of a proton with infinite 
mass. The only difference is the appearance of different masses in the various terms of hf,i and hpe- However, by 
appropriately redefining the length scales and dimensionless parameters involved in the problem, one can map /lee 
onto hf,i. This means in detail the replacement of Ag by Aee and of by in Eq. ( p2| , p3[ ) . Now we may use the 
analyticity of the virial coefhcient with respect to the interaction parameter. We may extend the result obtained for 
the electron- ion part Eq.(^) by analytical continuation to negative ^-values. Thus we have for the electron-electron 
contribution 

Bt = 2^\l,Q\U) + 2^AL||Q^(6e) . (56) 

The magnetic quantum virial function Q^{^ee) is given by Eq. js^) . Note that this series holds also for the ion-ion 
interaction if me is substituted by m^. However its contribution to the virial coefficient is negligible in the weak field 
limit. 

Let us briefly state the result for the exchange part of the electron-electron contribution. It may be obtained 



by introducing an additional factor (—1) in Eq.(B4 and 46) which takes into account the exact symmetry of the 



wavefunction. Then following the same steps as described in section HI we find 



^3 cosh(2a^^o^^ ^ a:^ cosh(2a;e) 

'cosh^(a;e) '^'^ 6 cosh^(xe) 



^ee ^ ~^-^ee ,2, (^ee) " 7^ Aee — . . 2/ i'^ee), (57) 



with the new magnetic exchange virial function 

'^''««") = £ foTiTJ) (ill" - ^"*'«' - ^' - ill" - ^"'>«"^ -")(¥)'■ <^«) 



The factor cosh(2a;e)/ cosh^(xe) in Eq.(|58|) is a result of the spins of the particles and can be calculated exactly. Again, 
one may check these results for the order fc = 0, 1 by comparison with the exact contributions given by Eqs.(p3 and 



C. Asymptotic properties of the virial function 



Let us now make an independent test of the above made statements. This investigation relies on two facts. First we 
consider the elctron-electron contribution only, then in the limit ^ ^ oo the quantum virial function (^) should be 
equal to the Wigner-Kirkwood expansion since ^ ~ Ti~^ . That means in this limit the plasma behaves essentially 
as a classical system. The second argument is that the electron-electron contribution may be obtained from the ion- 
electron contribution, and vice versa, by simple replacements of the interaction parameter as discussed in the previous 
section. Let us start by studying the higher order contributions fc > 6 to the magnetic virial function (truncated virial 
bmction), which read according to Eq.(^) as 

n'B(P\ /- V (fc ~ 3)C(fc - 2) + (fc - 5)C(fc - 4) / A ^ 

Q ^^) = ^^^ r(2 + fc/2) UJ ' ^ ^ 

with ^ < 0. The F- function can be represented by an inverse Laplace transform 
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1 f^+'°^ dt 



e 



With that and after rearranging the sum over k, Q'^(C) can be rewritten as 



(60) 



(?) 




r(4) \2) ^ r(9/2) V2 

In the following we make use of the relation 



;fflfiV-^ig2,Ci| . ,61, 



which gives then 



(0 



dt 



2nit^\2ViJ \ \2VtJ \ \2ViJ 



:C(2) /? 
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It useful to employ the asymptotic expansion of the tp function 

1 ™ B 

^(x) -Inx- — -^^-^+r™(x), (64) 

s— 1 

with the Bernoulli numbers i?2fc- Then we can perform the inverse Laplace transform and find the following asymptotic 
expansion of the truncated magnetic quantum virial function 

= (1 + 2C(3)) (IX-^JU ({2)) (lY-^Al + B,)^^'" 



r(9/2) ' ^'''V2/ r(4)V2 '7V2/ r(7/2) 



(i?2 + S4) U - ^^^iBi + i?6) 7 - -^^(^6 + Ss) 7 + o(r') • (65) 



2r(3)V2/ r(5/2) ' ' "'V2/ r(3/2)' " Vu r(i/2)' " °' 

Now we may conclude that the full magnetic virial function defined by 

5 



3 



(fc-3)C(fc-2) + (fc-5)C(fc-4) 

fe=2 

has the following asymptotic representation 



Remarkably, this procedure is accompanied by a term by term cancellation of the lower order contributions fc < 6 
coming from the Taylor expansion with those coming from the asymptotic expansion. The final expression may now 
be compared with the expansion , which was computed by Alastuey and Jancovici ]T^ . The linear term in the 
asymptotic expansion of Q^{£,) is the term proportional to ?i^-term of the Wigner Kirkwood expansion and coincides 
with that of Alastuey and Jancovici. In addition to that we have found higher order contributions proportional to 
and . With this derivation we have shown that that the magnetic virial function (53) has the correct asymptotic 
properties. This may be regarded as a strong support of the argument that the sum in Eq. (^9|) can be extended to 
the values of k=2,3,4, and 5, in order to obtain the desired result as given in Eq.([5^). 

Notice that from the Wigner-Kirkwood expansion follows the absence of the linear term in the Taylor expansion 
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Finally we give the asymptotic form of Q^{£,) for positive arguments. To establish this property, we first observe 
that the magnetic virial function obeys the following relation 



fc=2,4, 



(fc-3)C(fc-2) + (fc-5)C(fc-4 ) 

r(2 + k/2) 



(68) 



From this, it follows by using the representation of the C- function as an infinite sum and then carrying out the sum 
over k, that 



(69) 



where we have defined 

oo 

a^(0 = E2"Ml 

oo 

-E 



1 



(5 



7n^ 



5(2) 





r 


1 


n2 2! \2 










(1 


I 2/ ^2 


2! 





1 

3! 



(70) 



Now let us briefly summarize the properties of the magnetic quantum virial function. In Fig.l we have plotted Q^{£,) 
for both positive and negative arguments, i.e. for electron-ion and electron-electron interaction , respectively. It shows 
an asymmetric behavior. For opposite charged particles the magnetic quantum virial function increases exponentially 
at large ^, i.e. at low temperatures, due to the formation of bound states. While for like charged Q^{£,) increases 
linear at large ^. 

The behavior of the exchange magnetic virial function is shown in Fig. 2. In the quantum regime, at small ^, one finds 
a finite contribution to the thermodynamic funcitons. While E^{£^) decreases exponentially in the classical regime, 
i.e. at large ^-values. This result was also found in [T^ . 



IV. MAGNETIZATION AND MAGNETIC SUSCEPTIBILITY 



We now compute the magnetization in linear approximation (weak-field limit) and construct from this the mag- 
netic susceptibility. Thereby spin effects and orbital effects are treated on equal footing. Let us suppose that the 
magnetization may be divided into ideal and interaction contributions 



(71) 



We restrict ourselves to the magnetization of the electronic subsystem, since the magnetization of the subsystem of 
the heavy positive ions is negligible small. However the contribution of the electron-ion interaction is fully included 
into our calculation. The ideal magnetization M = —{l/Q,){dF/dB) may be calculated from 



M, 



id 



nkT-^ In 
oB 



2 tanh(a::e) 



(72) 



Evaluating this in the weak field limit, we get Landau's result for the sum of the spin magnetism and the diamagnetism, 
which reads 



Mr 



1 nh^e^PB 



6 



(73) 



The interaction part of the magnetization may be expressed in terms of the magnetic virial function. By taking the 
derivative of the full second virial coefficient with respect to the magnetic field one obtains 



Mi, 



7. 5 -j'^KQe + T 

6 to2 V4 " " 4 



B 

ee ^ ee 



37rnALi?°e 



(74) 



Here we have introduced a density expansion of the thermodynamic functions, that can be obtained from the fugacity 
expansion by an iteration procedure, as discussed in p^ , pO| . This expression may now be used to calculate the zero 
field magnetic susceptibihty (x = (O^Ml + Mint) / dB) ^Z^, with the result 
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X = XL (l + \n\lQ^{^,) + JnALQ^fee) - innXl^E^U) ~ \n\l^E^ {i,,)) . (75) 

The first term is Landau's result for the magnetic susceptibility, xl — (l/6)(7i?i^e^/?/TOg), of an ideal system in 
Boltzmann statistics, while the next terms describe the density effects. These effects contain the interaction of 
the particles as well as the deviation from the Boltzmann statistics. Fig. 3 shows the magnetic susceptibility as a 
function of the density parameter nA^g of the system for various temperatures. In Fig. 4 we have plotted the magnetic 
susceptibility as a function of the inverse temperature for various fixed densities. We find for fee < 1.2 a decrease 
and for ^^e > 1-2 an increase of the paramagnetic susceptibility. The transition from negative to positive corrections 
occurs at r ~ 2 X 10^ K. This non-monontonic dependence on the temperature is the result of two competiting effects. 
The first effect can be explained on the basis of an ideal quantum plasma. The exchange contribution of the order 
n^, which describes the first deviation from the Boltzmann statistics, decreases the magnetic susceptibility. On the 
other hand, the interaction between the particles tends to increase the magnetic susceptibility. This effect becomes 
dominant at low temperatures, while at high temperatures the exchange effects are dominant. 

We note that for f ^ 1, i.e. for T <C 2 x lO^K, the contribution from the positive interaction parameter (f > 0) may 
become very large due to its exponentiell increase with 1/T. The region where a considerable number of bound states 
are formed, requires a special treatment pc| ] . Clearly, this theory is restricted to the region in which \ x~Xl \ Ixl < 1 
is valid. 

Finally, we mention that the magnetization and magnetic susceptibility of an OCP can be derived from the results 
of the TCP ( ]7^ , [75| ). This limit is obtained by sending the mass of one species to infinity and the charge to zero while 
ensuring charge neutrality of the system. Then the magnetic susceptibility of an OCP in linear response reads 

X^'^'^ - XL (l + ^nXlQ^iU) - i^n\l^E\U) - ^nXl^E^ {U)) ■ (76) 

In the previous section, we have checked that this expression coincides with the Wigncr-Kirkwood expansion derived 
by Alastuey and Jancovici |p^ . 
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Figure Captions 

Fig. 1 Plot of the magnetic quantum virial function Q^{£,) ■ The positive branch > 0) corresponds to the 
electron-proton interaction and the negative branch < 0) to the electron-electron interaction. 
Fig. 2 Plot of the exchange magnetic quantum virial function E^{^). 

Fig. 3 Magnetic susceptibility as a function of the degeneracy parameter nX^^ for various temperatures (note that 

1^1 - (157000/r[X])'/'). 

Fig. 4 Magnetic susceptibility as a function of the of the coupling parameter F (inverse temperature) for various 
fixed densities. 
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